Quantum relaxation after a quench in systems with boundaries 
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We study the time-dependence of the magnetization profile, mi(t), of a large finite open quantum 
fsing chain after a quench. We observe a cyclic variation, in which starting with an exponentially 
decreasing period the local magnetization arrives to a quasi-stationary regime, which is followed 
by an exponentially fast reconstruction period. The non-thermal behavior observed at near-surface 
sites turns over to thermal behavior for bulk sites. Besides the standard time- and length-scales a 
non-standard time-scale is identified in the reconstruction period. 



Recent experimental progress in controlling ultracold 
atomic gases in optical lattices has opened new per- 
spectives in the physics of quantum systems. In these 
measurements the coupling in an interacting system can 
be tuned very rapidly, commonly denoted as "quench" , 
for instance by using the phenomenon of Feshbach reso- 
nance and the couplings to dissipative degrees of free- 
dom (such as phonons and electrons) are very weak. 
As a consequence one can study coherent time evolu- 
tion of isolated quantum systems. Among the fasci- 
nating new experiments we mention the collapse and 
revival of Bose-Einstein condensates [lj], quenches in a 
spinor condensate 0, realization of one-dimensional Bose 
systems Q and measurements of their non-equilibrium 
relaxation^. 

Concerning the theoretical side of quantum quenches 
here the first investigations had been performed on quan- 
tum XY and quantum Ising spin chains [H-Q before the 
experimental work has been started. The new experi- 
mental results in this field have triggered intensive and 
systematic theoretical researches, which are performed 
on different systems, such as ID Bose gases j8|, Luttinger 
liquids @ and others Besides studies on specific mod- 
els there are also field-theoretical investigations, in which 
relation with boundary critical phenomena and confor- 
mal field-theory are 

utilizeddlGl. 

One fundamental question of quantum quenches con- 
cerns the nature of the stationary state of this non- 
equilibrium quantum relaxation including the issue of 
thcrmalization and potential descriptions by Gibbs en- 
sembles. For non-intcgrable systems exact thermal- 
ization of stationary states was conjectured (l3j . how- 
ever the numerical results on specific systems are 
controversial 13l4l5l| . On the other hand integrable sys- 
tems are sensitive to the initial states and their stationary 
states are thermal-like being in a form of a generalized 
Gibbs ensemble @. 

Thcrmalization includes generically (i.e. away from 
critical points) an exponential decay of correlation func- 
tions in the stationary state on length and time scales 
that can be related to the correlation length and time of 
an equilibrium system at an effective temperature de- 



pending on the parameters of the quench 0, HH, G3- 
Some quantum systems do not thermalize completely and 
display a different behavior for correlation functions of lo- 
cal and for non-local operators, such that the former do 
not exhibit effective thermal behavior [l^]. An interesting 
issue not being addressed so far is the characterization of 
the non-stationary, that means not time-translation in- 
variant, quantum relaxation following a quench: Prepar- 
ing the quantum system in a non-eigenstate of its Hamil- 
tonian, how is thcrmalization achieved during the time- 
evolution? How do correlations develop in time towards 
the stationary (i.e. time translation invariant) state, is 
there a time dependent correlation length, etc.? 

Another important issue concerns quantum relaxation 
and potential thermalization in the presence of bound- 
aries. Theoretical studies of non-equilibrium quantum 
relaxation have focused on bulk sites up to now, but all 
real systems have a finite extent and they are bounded by 
surfaces and the physical properties in the surface region 
are considerably different from those in the bulk 18j . Ob- 
viously an interesting question is whether the time and 
length scales characterizing the stationary relaxation in 
the bulk is altered in the vicinity of the boundary, and 
how thermalization is achieved there. 

In this paper we will address these two issues: The 
non-stationary quantum relaxation after a quench and 
the effect of boundaries. For this we focus on a compu- 
tationally tractable model for a quantum spin chain and 
study the relaxation of profiles of observables in the early 
time steps as well as their behavior in the long-time limit. 
We also address the behavior in large, but finite systems 
and study the consequences of the recurrence theorem. 

The system we consider in this paper is the quantum 
Ising chain defined by the Hamiltonian: 
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in terms of the Pauli-matrices af' z 



at site I. In the 
non-equilibrium process the strength of the transverse 
field is suddenly changed from h (t < 0) to h (t > 0). 
The Hamiltonian in Eq.((T]) can be expressed in terms 
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TABLE I: Decay exponent of the off-diagonal (longitudinal) 
magnetization in the initial (equilibrium) period. 





ho = h c 


ho > h c . 


bulk 
boundary 


1/8 
1/2 


1/2 
3/2 



of free fermions[19|, which is used in studies of its non- 
equilibrium properties 0, [l6|. The bulk transverse mag- 
netization, which is a local operator, has non-thermal 
behavior [a fill flTT] . whereas the bulk (longitudinal) mag- 
netization, erf, which is a non-local operator, has effective 



thermal behavior 
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Here we concentrate on the latter 
quantity and study the time-dependence of its profile, 

mi{t) = lim 6 ^ 0+ b (H 0) \<rf(t)\*i 0) )b, wh ere |^ 0) ) h is the 
ground state of the initial Hamiltonian (TTJ) in the pres- 
ence of an external longitudinal field b. According to [20[ 
this can be written as the off-diagonal matrix-element of 
the Hamiltonian ([T]): 



(2) 



Here 1^1°^) is the first excited state (which is the ground 
state of the sector with odd number of fermions) of the 
initial Hamiltonian (t < 0). In the ordered phase, ho < 
h c = 1 where mi(t < 0) > 0, ) is asymptotically 

degenerate with the ground state, \Wq °'). For ho > h c 
the magnetization vanishes as vrii (t < 0) ~ L~ x with the 
system size for t < 0. The decay exponent, x, is different 
at the critical point, h = h c , and in the paramagnetic 
phase, h > h c , as well in the bulk (l/L = 0(1)) and at 
the boundary (l/L — >• 0), see Table [U 

To calculate the magnetization profile in Eq.([2t we 
have used standard free-fermionic techniques fl9l 1 2 lj . For 
the surface site, 1 = 1, most of the calculations are an- 
alytical, whereas for / > 1 numerical calculations have 
been made for large finite systems up to L = 384. 

We have performed quenches for various pairs of trans- 
verse fields, ho and h and calculated the time-dependence 
of the local magnetization at different sites, I < L/2. The 
results depend primarily on whether the system before 
and after the quench is in the ordered (O) or disordered 
(D) phase, see Fig Q] for different combinations of O and 
D. One can identify different time regimes that can be 
interpreted in terms of quasi-particlcs, which arc emitted 
at t = 0, travel with a constant speed, v = v(h, ho), and 
are reflected at the boundaries. 

As argued in Ref.[ll[ orn y those quasi-particles are 
quantum entangled that originate from nearby regions 
in space, others are incoherent. When the latter arrive 
at a reference point I they cause relaxation of local ob- 
servables (such as magnetization). Here we extend this 
picture by noting that in a system with boundaries the 
same quasi-particle can reach the point I twice (or more) 
at different times after reflections. This induces quan- 




FIG. 1: (Color online) Relaxation of the local magnetization, 
log mi(t), at different positions in a L = 256 chain after a 
quench with parameters: a) ho = 0.0 and h — 0.5 (O —5- O) 
b) h = 0.5 and h = 1.5 (O -> D) c) h = 1.5 and h = 0.5 
(D ->■ 0)d) h = 1.5 and h = 2.0 (D -J- D). 



turn correlations in time signalized by the reconstruction 
of the value of the local observable. In the following we 
analyze the different regimes of the relaxation. 

In the free relaxation regime: t < ti = l/v only inco- 
herent quasi-particles pass the reference point resulting 
in an exponential decay of the magnetization (cf. Fig. [IJ : 



mi(t) = m(t) « A(t) exp(-t/r), t < ti 



(3) 



with an oscillating prefactor, A{i). In the regime h > h c 
and ho < h c we have A{t) ~ cos(at + b), thus m(t) 
changes sign. On the other hand in the other parts of 
the phase diagram m(t) is always positive, i.e. A(t) ~ 
[cos(at + b) + c], with c > 1. The characteristic time- 
scale, r = r(h,ho), is the relaxation or phase coherence 
time, which is extracted from the numerical data. The 
exponential form of the decay in Eq.Q indicates ther- 
malization, at least for bulk sites, which is in agreement 
with the similar decay of the autocorrelation function. 

In the quasi-stationary regime: ti < t < T — ti , T = 
L/v, two types of quasi-particles reach the reference point 
I: type 1 passed I only once at a time t' < t and type 
2 passed it twice at two times t' < t" < t with a reflec- 
tion at the nearby boundary between t' and t" . These 
two types interfere, resulting in a comparatively slow re- 
laxation (cf. Fig. [1} . Deep inside the ordered phase the 
quasi-particles can be identified with kinks moving with 
a speed ±v [22[ and in the regime (; < ( « T half of the 
quasi-particles reaching the site / are of type 1 (flipping 
the spin at / once) and half of them type 2 (flipping it 
twice), leading to a quasi-stationary relaxation. 

The magnetization profiles for fixed times t < T/2 are 
shown in Fig|2] for the same quenches as in FigJTJ For 
sufficiently large I the quasi-stationary magnetization has 
an exponential dependence, such that comparing its value 
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TABLE II: Correction to the quasi-stationary behavior for 
the surface magnetization in different domains of the quench. 





ho < h c 


h Q > h c 


h < ho 
h > h 


t^ 1 cos(at + b) 
t~ 3/2 cos(at + b) 


L~ 3/2 [cos(at + b) + c], c>l 
t-^ 2 [cos(at + b) + cL^ 2 } 



at two sites, l\ and I2, we have 

m hih) / m h {t 2 ) ~ cxp - l 2 )/£,} 



(4) 



with oscillating prefactors. 

In the limits L — -> 00 and t — > 00 one can define a 
quasi-stationary limiting value which will be denoted by, 
to/. For the surface site we have the exact result 



_ _ (l_ fe 2 )(1 _^ ) l/2 

f — hh,Q 



h 0l h < I , 



(5) 



and zero otherwise. Note that the non-equilibrium sur- 
face magnetization has different type of singularities for 
h ->• 1~ (h < 1) and for ft ->• 1~ (h < 1). We have 
analyzed the correction term, A(t,L) = mi(t) — mi, and 
its asymptotic behavior is summarized in Table [TT] in the 
different domains of ft and fto- These corrections are in 
power-law form, which signals that the relaxation of the 
surface magnetization has non-thermal behavior. 

For ( > 1 we observe that to; is monotonously de- 
creasing with / and thus m; > for ho , ft < 1 and zero 
otherwise. The correction terms are identical with those 
given in Table [TT] so that a finite distance, I, the local 
magnetization has non-thermal behavior. 

In the reconstruction regime: T — ti < t < T more and 
more quasi-particlcs of type 2 reach the reference point, 
which implies, within a kink-picture, that incoherent 
spin flips in the past are progressively reversed by quasi- 
particles returning to the site / after reflection. For mono- 
disperse quasi-particles (velocity v) one would expect a 
T-periodicity and thus to; (t) = to; (T — t) , i.e an expo- 
nential increase in t with a growth rate similar to the 
initial decay rate. Indeed we find 



to; (t) = m(t) w B(t) exp(i/i 



T-ti<t<T, (6) 



which is practically position independent and where the 
growth rate of r'(h, h ) depends on the conditions of the 
quench, being approximately proportional to r(h, ho): 
t/t' = 0.883 ± 0.002. It turned out to be useful to 
measure the cross-over time, t = T/2, which is defined 
as the crossing point of the two asymptotic regimes: 
Aexp(— t/r) = Sexp(t/T'), where A and B are aver- 
aged prefactors. During the cross-over time the quasi- 
particles travel a distance, L/2, thus their speed is given 
by: v{h, ho) = L/2t, which can be measured accurately. 
We have noticed, that for h < 1 the speed is proportional 
to h: v(h,ho) = ha(h,ho), where a{h,ho) is practically 
independent of ho and has just a very week dependence 
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FIG. 2: (Color online) Non-equilibrium magnetization pro- 
files, logm;(t), at different times after a quench with param- 
eters given in Fig[T]for L — 384. From the asymptotic values 
of the slopes one can measure the correlation length. 



on h close to ft. = 1. The typical values are in the range 
a(h,ho) ~ 0.86 — 0.88. For h > 1 the speed is practically 
constant and has no h dependence. 

Approximate periodicity with T starts for t > T, when 
quasi-particles start to be reflected second time and the 
spin-configuration of the system becomes approximately 
equivalent to that at t — T '. 

The time- and length scale, as defined in Eq.© and 
Eq.Q, respectively, as well as the characteristic quasi- 
particle speed v(h, ho) = £/t, can be extracted with high 
numerical accuracy from our data for the magnetization 
profiles, typically with a precision of 3 — 4 digits. Com- 
plementary calculations of the autocorrelation function 

Gi(t) = (4 0) 



relation function, Ct(r) 



r(°)\ 

<*o (0) l 



i {t)crf (0)\Wq and the equal-time cor- 



l+r 



(t)af(t)\4 0) ) show 



that they yield the same correlation time and length, but 
with less accuracy. Based on our results for the profiles 
we have conjectured possibly exact results about the re- 
laxation time, as discussed below. 

The relaxation time r(h, h ) is divergent at two 
points: i) at the stationary point, h = ho, where 
r(h, ho) ~ (h—ho)~ 2 and ii) for small h, where r(h, ho) ~ 
h , which can be derived perturbatively. For ho = the 
two singularities merge at h = 0: r(h, ho — 0) ~ h~ 3 . 

To obtain information about r(h, ho) away from the 
singularities we consider a quench from the fully ordered 
initial state (ho — 0) first. A quench into the disor- 
dered phase (h > I) yield to high numerical accuracy 
r(h > l,ho = 0) = 7r/2, i.e. independent of h. For 
a quench into the ordered phase (h < 1) we introduce 
f(h, ho = 0) = h 3 r(h 7 h = 0) to get rid of the singu- 
larity at h = 0. In the limit h — > we obtain f(h = 
0, ho = 0) = 37r/2, and for ft, > we consider the ratio: 
y T (h) = Af(/i)/Af(0) with Af(h) = f(h) - f(l) and 
compare it with a similar expression for the correlation 
length yt(h) = A|(ft)/A/|(0) with Af(ft) - |(ft) - |(1), 
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where £(ft) = £(ft)ft 2 . The two ratios y T {h) and jr(ft), 
as shown in Figj3£i, are almost indistinguishable. Since 
£(ft) = -l/log((l + \/l - /i 2 )/2) is known exactly 0, the 
relaxation time for a quench from an ordered initial state 
(fto = 0) can therefore be estimated very accurately, if 
not exactly, by the relation y T (h) = 2/^ (ft). 




FIG. 3: (Color online) Left: The ratios y* = Af(ft)/Af(0) 
and y T = Af (ft)/Af(0) for a quench from Tio = as a func- 
tion of h. The curve y ? derives from the exactly known 
form for £(/i), see text. Right: The ratios y T (h,ho) = 
Af (ft, fto)/Af (0, 0) = r(ft, ft ) • ft(ft - ho) 2 /V - (1 - h f/2 
for various fto as a function of ft. 

Starting from a partially ordered initial state (0 < fto < 
1) we define f(h, ho) = h(h — fto) 2 r(ft, ho) and find to 
high numerical accuracy that the limiting value at h = 1 
is given by: f(h = 1, ho) = 7r(l — ho)/2. Away from h = 1 
we study the ratio y T (h,ho) = Af (ft, fto)/Af (0, 0) with 
Af(ft, fto) = f (ft, fto) — f (1, fto) which is identical to y T (h) 
for fto = and which is plotted in Figj3]3 for different 
values of fto . The curves for all values of fto are quite close 
to each other, and at ft = 1 they all have a singularity, ~ 
vl — ft. Therefore one obtains a very good estimate for 
the relaxation time from f(ft, fto) by y 7 (ft, fto) ps y T {h) = 
y£ (ft), which is given in an analytical form (see above). 

The thermal-like stationary state can be character- 
ized by an effective temperature T e // [HI which is de- 
fined through the condition, that the relaxation time in 
the stationary state after a quench, r(ft, fto), and the 
equilibrium correlation time at temperature T = T e //, 
TT(h,T), are identical. Using the analytic result at the 
critical pointjlH: ry(ft = 1,T) = 8/(wT) we arrive at 
T e ff(ho,h = 1) = 16(1 — fto)/7r 2 , which is compati- 
ble with the numerical data in Ref.[3|. In the ferro- 
magnetic phase, ft < 1, and in the limit T -C A (ft) , 
A(ft) being the gap, the relaxation time is given by [221 ] : 
r T (ft < 1,T) ps (2/(7rT))e A / T , which for |ft - ft | < 1 
leads to: T eff ss -A(ft)/(21n|ft - ft |). 

To summarize we have identified different regimes in 
the non-equilibrium relaxation of the magnetization pro- 
files of the quantum Ising chain with boundaries, which 
can be explained in terms of quasi-particles that are re- 
flected at the surfaces. For sites at or near the sur- 
face non-thermal behavior is observed, manifested by a 
power-low relaxation form. For bulk sites a cross-over 
to thermal behavior is found, with exponentially decay- 
ing correlations, defining a relaxation time and a correla- 
tion length that is identical in semi-infinite and in infinite 



systems and which obey presumably exact relations con- 
jectured on the basis of the numerical data. In a finite 
system an exponentially fast reconstruction of the local 
magnetization is observed, involving a time-scale, r', and 
characterizing an approximately periodic dynamics. 

Several results for observables displaying thermal be- 
havior in the bulk are expected to be valid also in other, 
even non-intcgrablc spin chains: Absence of thcrmaliza- 
tion at the boundaries, identity of correlation time and 
length in infinite and semi-infinite systems and an expo- 
nentially fast reconstruction in finite systems. 
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